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We consider an interacting pion gas in a stage of the system evolution where thermal but not 
chemical equilibrium has been reached, i.e., for temperatures between thermal and chemical freeze- 
out Tther < T < T c hem reached in Relativistic Heavy ion Collisions. Approximate particle number 
conservation is implemented by a nonvanishing pion number chemical potential [i w within a di- 
agrammatic thermal field theory approach, valid in principle for any bosonic field theory in this 
regime. The resulting Feynman rules are derived here and applied within the context of Chiral 
Perturbation Theory to discuss thermodynamical quantities of interest for the pion gas such as the 
free energy, the quark condensate and thermal self-energy. In particular, we derive the fi^ 7^ 
generalization of Luscher and Gell-Mann-Oakes-Renner type relations. We pay special attention to 
the comparison with the conventional kinetic theory approach in the dilute regime, which allows 
for a check of consistency of our approach. Several phenomenological applications are discussed, 
concerning chiral symmetry restoration, freeze-out conditions and Bose-Einstein pion condensation. 

PACS numbers: ll.10.Wx, 12.39.Fe, 25.75.-q. 



I. INTRODUCTION 

One of the ongoing research lines in Heavy Ion Physics is the thermal and chemical evolution of the expanding 
hadronic gas. Roughly speaking, the accepted picture is that the evolution of the cooling system reaches chemical 
freeze-out before the thermal one, so that when hadrons fully decouple the chemical potentials associated to particle 
number conservation are not zero. The chemical composition of the gas can be determined experimentally by looking 
at the relative abundances of the different hadron species and their spectra [IJ [H [31 El [5] . The presence of such a 
chemically not-equilibrated phase is more likely to exist for higher collision energies such as those in RHIC or LHC 
than for SPS or AGS experiments [4 . For the pion component, different estimates based on local thermal equilibrium 
and particle spectra analyses predict fi^ ~ 50-100 MeV at a thermal freeze-out temperature T t her ~ 100-120 MeV, 
with chemical freeze-out taking place at about T c h em ~ 180 MeV [H [3J 03 GO [7J [B] . On the other hand, the plasma 
is almost electrically neutral, so that it is a good approximation to keep vanishing charge or isospin chemical pion 
chemical potentials. 

For low and moderate temperatures, the dominant component is the pionic one. In that phase, the mean free path 
of pions is small compared to the system size, so that local thermal equilibrium prevails [2J [5J [ID]. On the other 
hand, the chemical relaxation rate through inelastic -kit t» nnn-K processes is very small [U [TT] due to a strong phase 
space suppression. Therefore, in the range of temperatures T t her < T < T c h em <>T c , with T c the chiral restoration 
critical temperature, the system is in thermal equilibrium and dominated by elastic collisions so that fi v =/= 0. In 
that temperature range, it is valid to use the theoretical framework of Chiral Perturbation Theory (ChPT) and it 
is also reasonable to adopt a dilute gas description, since the mean particle density is small. In addition, neglecting 
dissipative effects such as viscosities, entropy is conserved in the evolution. 

The system described above, i.e., a pion gas with ^ is the one we will consider here. Clearly, it is an 
oversimplified version of the real hadron gas, but we will take it as a physically relevant working example for our present 
analysis. So far, pion number chemical potential effects in such a system have been incorporated basically in two ways. 
One of them is the limit of free particles (where one has actually exact particle conservation) used for the evaluation 
of the partition function, including resonances explicitly [2,8. This allows, via entropy conservation requirements, 
to describe rather accurately the isentropic dependence ^{T) in the range of temperatures of phenomenological 
relevance indicated above. The other one is to use kinetic theory arguments to include the fi^ 7^ dependence 
directly in the distribution function. The latter has been followed for instance in the calculation of the thermal 
width 9 , in the evaluation of transport coefficients [TH [T3] or in the virial approach for low densities [2] . Finally, 
it is worth mentioning that there are phenomenological analyses, like that in |15j for the dilepton rate, where the 
same prescription is followed, i.e., replacing the distribution function, but for propagators at the diagrammatic level, 
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inspired on the nonequilibrium formulation of thermal field theory [16 . 

We will be interested in a diagrammatic formulation of this system, i.e, we will derive the Feynman rules to be 
used when approximate particle number conservation is valid. The Feynman rules of thermal field theory with exact 
conserved charges can be obtained straightforwardly [T7] but this is a completely different situation, since particle 
number is not exactly conserved in an interacting bosonic field theory (is only conserved in the free case) and therefore 
there is not a local number operator to be added to the lagrangian in the usual way. This will also be related to the 
impossibility to define a proper Matsubara imaginary-time formalism. The motivation for our field-theory description 
is twofold: first, it will provide a formal proof of the consistency and validity of the different prescriptions used in 
the literature and mentioned in the previous paragraph. Second, it will allow to deal in a natural way with pion 
interactions when fj, w ^ 0, which is particularly interesting in order to describe corrections to dynamical quantities 
such as the thermal pion self-energy, but also to evaluate the effect of interactions in thermodynamical observablcs. 

The paper is organized as follows: in the first part we will describe our formalism, based on holomorphic path 
integrals, which naturally leads to the relevant Feynman rules. Although the results in that part are actually valid 
for any real scalar field theory provided one neglects the contributions of number-changing scattering processes or, in 
other words, if the gas is dilute enough, we will be primarily interested in the pion gas, where chemical nonequilibrium 
is actually reached during the expansion. We explain more clearly this physical motivation in section |TT| w here we 
discuss the relevant distribution function to describe the system. The Feynman rules we derive (section III) had not 
been considered before in the interacting case and, as we will see, they are really meaningful only in the real-time 
formalism of thermal field theory. The second part (section |IV| ) deals with the application of our formalism to the 
pion gas. We will analyze corrections both in thermodynamical (free energy, entropy, particle number and quark 
condensate) and dynamical (thermal mass and width) observablcs, the former being understood as a generalization of 
the usual thermodynamical variables during the chemical nonequilibrium phase. We also compare to previous works 
in the literature and we discuss several phcnomcnological consequences regarding chiral symmetry restoration, Bose- 
Einstein condensation of neutral and charged pions, as well as thermal and chemical freeze-out. Appendices |A"| and [B] 
contain detailed results used in the main text about holomorphic path integrals and thermal propagators respectively. 
In particular, in Appendix [B] we discuss some relevant differences between the case of particle number chemical 
potential considered here and the more usual one associated to the electric charge exact conservation, concerning 
especially the way in which KMS boundary conditions are broken. 



II. PHYSICAL MOTIVATION 



As stated in the introduction, we are interested in describing the physical system constituted by a pion gas in 
expansion, during the time when the number of pions is approximately conserved. This is the case of the pionic 
component of the hadronic gas produced after a relativistic heavy-ion collision. The pionic component is the dominant 
one in the hadron gas around thermal freeze out pQ , although considering additional degrees of freedom in the gas 
(kaons, etas, nucleons and resonances) and the interactions among them is relevant at temperatures close to the chiral 
phase transition. We will not consider those extra components here, although their inclusion in the chiral framework, 
together with the corresponding additional chemical potentials (strangeness and baryon number) is a feasible extension 
of this work. Unlike other treatments 0, where it was shown that it is a reasonable approximation to introduce in 
the partition function all the states (asymptotic states as well as resonances) up to a given energy as free degrees 
of freedom, in our approach the resonances present in the pion gas, the /o(600)/ct and the p(770), are generated 
dynamically by means of unitarization methods so the actual degrees of freedom in the lagrangian are only pions. It 
is however important to mention that, even when introduced as explicit degrees of freedom, the processes p ^ tttt 
and a rnr do not restore chemical equilibrium in the pionic component, because /i p = )i a = the truly relevant 
particle-changing process being tttt ^ tttttttt. When pions and resonances are in chemical equilibrium with respect 
to each other we talk about a relative chemical equilibrium, since it is possible to choose their chemical potentials to 
maintain it, whereas absolute chemical equilibrium is only possible for fi n = [S]. 

The evolution of the pionic fireball can be divided into three stages as it cools down according to the corresponding 
temperature ranges [2J: (I) Tc^em <T <T C1 the pion gas is produced after hadronization from a quark-gluon plasma 
phase and it is in full statistical equilibrium (thermal and chemical). (II) T t her <T< T c / lem , the mean- free path of 
elastic collisions, A e ; is smaller than the typical size of the fireball, R ~ 5 — 10 fm [18], so that thermal equilibrium 
is maintained, whereas the mean-free path of inelastic collisions \i n is larger than R so the total number of pions 
N-n = N n o + N n + + N„- remains approximately constant 1 and a finite chemical potential associated to N v builds up, 
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so that the system is out of chemical equilibrium. (Ill) T < T t her, A e / is larger than R so the pions stop interacting 
and their momentum distributions become frozen, so that thermal equilibrium is lost. 

We shall analyze the phase II of the evolution, where the total number of pions is approximately conserved (and the 
gas remains dilute enough) so that the introduction of a finite pion chemical potential \x v associated to N T is necessary. 
The existence of such a chemically not equilibrated phase during the fireball space-time evolution is supported by 
several phenomenological results. For instance when analyzing experimental data from the NA44 collaboration, it 
was shown in [3] that in order to fit properly the pion spectrum at low transverse momentum in PbPb reactions, one 
needs to introduce a finite chemical potential of order \i v ~ 60 — 80 MeV at thermal freeze-out. A similar conclusion 
is reached in [5] for RHIC AuAu collisions. In addition, the analysis of total particle yields and yield ratios for SPS 
and RHIC energies are fitted with values of fugacities compatible with the pionic component being significatively out 
of chemical equilibrium [Bj |7J . 

Neglecting electromagnetic interactions, the pions are described by neutral scalar fields. For a neutral boson field 
theory, particle number is conserved only in the free case. Our aim here is to provide a field-theory description of the 
non-equilibrium state corresponding to phase II. In this respect, it is important to remark that there are fundamental 
differences between total particle number and charges which are exactly conserved by the dynamics, such as the 
net electric charge or baryon number in QCD. In a pion gas, the total number of pions is expressed in terms of 
individual number operators as N — N„o + N v + + N n - whereas electric charge, or the third isospin component 
(they are equivalent in the pion SU(2) case) is measured by N n + — N„- [HI [20]. The main difference is of course 
that charge is exactly conserved in the second case, which implies several important consequences: first, from the 
field-theoretical point of view, in the charge case there is a local charge operator Q written in terms of the field and 
its derivatives, which allows for a straightforward derivation of the corresponding Feynman rules, adding the usual 
HqQ term to the lagrangian |17l 121] . However, that is not the case for particle number, which instead has a natural 
formulation in terms of canonical creation and annihilation operators. That is the main reason why we will develop 



our holomorphic representation of the path integral in section III Second, in the fj^ case we are really facing a 
nonequilibrium description, which is only consistent if fi v and T are not independent parameters, the function /^(T) 
parameterizing the deviations from chemical equilibrium and vanishing at T — T c h em . This signals the end of phase 
II, or its beginning if we think in terms of proper time, as inverse of temperature evolution in a hydrodynamical 
description (remember that in phase II local thermal equilibrium is assumed). The form of Htt(T) has to be fixed by 



additional physical assumptions. We will rely here (see section IV B I on the isentropic condition stating that the ratio 
of entropy density to pion density s/n remains constant along the chemical evolution, which has phenomenological 
support [H [3] • Finally, we remark that these differences between the charge \iq ^ case and the pion number /x w ^ 
one translate into a different way in which the KMS boundary conditions characteristic of equilibrium are broken. 
We discuss this issue in detail in Appendix 15] Throughout this work we will take hq — 0, which corresponds to an 
electrically neutral pion gas, which seems to De well supported by the phenomenological values of the fugacities [TJ. 

Since the system in phase II is in thermal equilibrium and there is an approximate conserved operator N with 
chemical potential fi v associated, the appropriate non-equilibrium partition function is 

Z p {t) = It , (1) 

where quantities with a tilde will refer always to the nonequilibrium /i,^0 case throughout this paper. By including 
source terms we can then derive thermal correlation functions. Note that Zp is independent of the position in space 
(we consider an homogeneous system), but it actually depends on (proper) time t during the gas expansion through 
temperature (3{t) = l/T(t) and the chemical potential fi(T(t)), along the lines discussed in the previous paragraph. 
The validity of the out-of- (chemical) equilibrium distribution function in will be subject to times t < tu, where 
tu is the duration of phase II. It is in this nonequilibrium effectively time-dependent situation that our results for the 



partition function and thermodynamical observables (see section [IV A I have to be understood 



For t < tu inelastic processes are scarce at temperatures T t h er < T < T c ^ em , therefore if at some time t\ the system 
is in a state with well defined number of particles equal to N, N\n(ti)) — N\n(ti)), then at another time t 2 with 
£2 — ti < tu, iV| 77,(^2) ) — N\n(t2)), and thus in Heisenberg's picture N(ti) ~ Nfa) and from Heisenberg's equation 
idN(t)/dt = [N(t),H] we infer: 

~ i(N(t 2 ) - N(h)) = I [N(t),H]dt => [N(t) , H) ~ , for t x < t < t n , (2) 



7T7r fc? 7T7r7r7r is Ti n ~ 200 fm 0, the typical mean velocities at those temperatures being v ~ c. See also our comments in section IV C 
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In the following section, the condition [N(t),H] ~ (valid for times in phase II) will be used to derive a field-theory 
description of the system based on the holomorphic path-integral representation of the partition function provided by 
Q. 

III. FORMALISM: CHEMICAL POTENTIALS FOR NEUTRAL BOSONS 

As we mentioned in the introduction, the operator N has a non-local representation in terms of the field operator 
|22j , so that the appropriate representation for this operator is instead in terms of creation and annihilation operators. 
The holomorphic path-integral representation uses these convenient operators, and its main ideas can be found in [23 . 
We will give the essential steps of the derivation for our system in this section, more technical aspects being relegated 
to Appendix [X] As we saw in the previous section, the generating functional of thermal correlation functions will be 
constructed from the non- (chemical) equilibrium partition function ([lj. In order to simplify further the discussion, 
we will consider first a quantum- mechanical gas of Bose particles and next we will extend it straightforwardly to the 
QFT case. Throughout this section and for simplicity we use the notation /j, for the chemical potential associated to 



particle number (/v in the pion gas case, which we will analyze extensively in section IV I. 

Let us consider then a single- frequency quantum oscillator (free Hamiltonian) coupled to an external force j(t). 
The Hamiltonian and number operators are then: 

H = l ~p 2 + -0J 2 f - j(t)q = H - j(t)q 

= |(ata + afit)--^=(at + 6)j(t) (3) 

N = a f a (4) 

where q and p are respectively the position and conjugate momentum operators (whose role will be played by the 
field and its conjugate momentum) and the creation and annihilation operators arc defined in the usual way: 

a = —i= (p — iuoq) a) = -={p + iuq) (5) 



2uj V2w 

satisfying canonical commutation relations [a, a 1 *] = 1. 

In the holomorphic representation |23] traces of operators are evaluated in the space of complex analytic functions 
of one complex variable z and creation and annihilation operators act on this space as: 

cl^ i i z , a i— > (6) 
oz 

We have included in Appendix [X] some of the technical details to perfor m th e relevant calculations in this formalism. 



In particular, the partition function for any Hamiltonian H reads, from (A7) 



Z,= j e -** (z\e-^»-^\z) (7) 



Now, if the number operator is approximately conserved, then [H, N] ~ and equation ([7|) can be recast, by 
inserting the identity once, as: 

ziri J 2ni 

This is the key step of the derivation, since it contains our main approximation, which is equivalent to consider only 
up to two-particle states in the trace Q. Therefore, it is physically appropriate to describe a dilute regime where 
elastic collisions dominate and particle number is approximately conserved. 



Now, the first matrix element in (|8| can be calculated directly, using (A14l with j = 0, i/ = i,- — i/3 and uj = —fx: 
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(z\e 



exp (/ 



(9) 



so that, using (A4) we arrive to: 



Z B = 



dz dz 
2iri 



(ze^\e-^ H 



(10) 



From this representation of the partition function we define the corresponding generating functional (in the QM 
case): 



zM 



dz dz 



(zeMe-^ 11 -^ 



2) 



(11) 



so that correlators of any function of the position operator q (the field operator in the QFT case) can be expressed in 
terms of functional derivatives of Zp [j] with respect to j at j = in the usual way. 

We will now proceed to the evaluation of Zp\j] when the Hamiltonian is the free one plus the source term, i.e, 
H = Ho — jq in Q3J) - Then, as usual, by functional derivation we will get the generating functional for the interacting 
case. We first separate the normal-ordered part as customarily, i.e, Hq = co/2 + Loa)a where the first term is the 
vacuum energy. Therefore, we have: 



dz dz 
2ni 



U (ze^,z;-if3) 



(12) 



The function W is defined in (A10) and for the present case, its expression is given in (Al4l-(A15l. For its evaluation, 
we have considered, as detailed in Appendix [X] the complex time contour showed in Fig |1| joining the points ti and 
ti — i(3 with a € [0, /?], which contains the usual real-time and imaginary-time paths of Thermal Field Theory and 
satisfies the usual requirements for the path integral to be well defined, i.e., Imt is monotonically decreasing along 
the contour |17j . The imaginary-time contour runs in a straight line from tj = down to — i(3 and is denoted as C4, 
while the C\ and C'2 are t h e path s used in the real-time formulation (see below). 



Now, replacing in (A14l-(A15l z 



ze 



z, t 



f 



ti — i/3 is equivalent to replace: 



/3->/3 = 



0(1 



(13) 



except for the (3 appearing in the C contour. With this replacement one can follow the same steps as in [SS] for the 
evaluation of the remaining integral in (12 1. Namely, one goes back to the discretized version of the path-integral (see 



Appendix |A| , uses again ( A13 1 with the modified A matrix and finally takes again the continuum limit. The final 
result is: 



t\ — ia 



13 



Im t 



-ti 



Ret 



FIG. 1: Complex time contour including real and imaginary time paths, used in the derivation of the /1 /; Feynman rules, 
where a £ [0,/3]. 
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Z$\j\ = Z»cxp | -- / dt dt' j(t)G F (t - t')0) 



(14) 



with the free partition function 



and the free propagator: 



G(t) 



p -/3w/2 

*0 i_ e -/9(w-^) 



l(l + n(w- ( u))+c^ |t| n(cj- /x) 



(15) 



(16) 



where the Bose- Einstein function: 



n(:r) 



(17) 



so that n(uj — /1) is the free distribution function at nonzero /1 in (B7l for a particle of positive energy. Note that 
we must restrict to /i < lo in order that the previous expressions for the partition function are well defined (sec also 
comments in Appendix m» . The upper limit would correspond to Bose-Einstein condensation (see below) . In the 
above propagator, \t — i'|has to be understood in terms of the relative position of times t and t' with the path routing 
shown in Fig[l] 

The result ( 14 I for the quantum-mechanical case for fj, 7^ is one of our main results. Its importance relies on the 
fact that we can now easily construct the generating functional in the interacting case, say H = Hq + V(q) with V 
the potential, in the usual way, i.e., by expanding formally in series of V and writing every term in the expansion in 
terms of functional derivatives of Zp[j] with respect to j. From there, the extension to a QFT for a real scalar field 2 
with lagrangian density: 



C 



-0 2 -V(0)-j0 



(18) 



is given by: 



Z p \j] = Z° p exp -i / d 4 xV 



iSj(x) 



exp 



d 4 x 



d 4 x'j(x)G(x - x')j(x') 



(19) 



where f c d 4 x = f„ dt J d 3 x. 



The generating functional (19) for the interacting case at [i 7^ and the corresponding Feynman rules which we 



discuss below constitute central results of this paper and, to the best of our knowledge, they had not been considered 
before. It is valid for any scalar theory, provided one works in the regim e wh ere clastic collisions domi nat e and particle 
number is approximately conserved. The propagator G appearing in (19) is the generalization of (16) to the QFT 
case when uj — > E p , the particle energy, and therefore we will restrict in the following to /1 < m. Recall that the 



QFT generalization of (3 in (|13[) is (3 V in (B13l. The explicit expression of the propagator coincides, as it should, with 



the free two-point function ( |B8| -(B11 1 derived in Appendix |B| directly within the canonical formalism at /1 7^ for 
t € K.. In this sense, one could somehow expect that the generalization to /.i 7^ of the generating functional is the one 
given in (19 1, although there was no rigorous proof available in the literature. We insist that the usual field-theory 



derivation for the case of an exactly conserved charge is not applicable here. 

Next, we will discuss the Feynman rules needed for diagrammatic calculations. The /i =/= case for approximate 
particle conservation is essentially a nonequilibrium situation, as commented several times before and thus it presents 
many subtleties to be borne in mind. One of them is the impossibility of defining properly a Matsubara or imaginary- 
time formalism, which is related to the way in which the KMS conditions are broken. We will separate this discussion 
from the real-time case, where a suitable formulation is possible, at least to the order we are considering here. 



2 We remark that in the field-theory case, whereas the Hamiltonian can be expressed as a space integral of a local field operator, that is 
not the case for the number operator when infinite frequencies appear. This is only possible for exactly conserved currents. 
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A. Imaginary-time formalism 



The Imaginary Time Formalism (ITF) corresponds to the choice ti = a = for the contour in Fig T] so that one 
is left only with purely imaginary times t = —it with r € [0, 0\. At /i — 0, this formalism is usually best suited to 
deal with thermodynamic quantities such as the free energy, while retarded Green functions can be derived from it 
by analytic continuation in the external frequencies [M]. However, the \i ^ propagator shows a distinctive feature 
that complicates diagrammatic calculations, generating in some cases ill-defined results. The origin of the problem is 
the way in which the standard equilibrium KMS periodicity conditions are broken. As explained in Appendix [B] to 
which we refer for notation, in our case (particle number conservation valid both for neutral and charged bosons) we 
have, in the mixed representation of the propagator, A^(r + (3 p ,p) = A^(r, p) instead of the familiar KMS condition 
Ar(r + P,p) — At(t,p). This momentum-dependent periodic condition makes it impossible to d efine properly a 
Matsubara representation in Fourier space, which can only be done for r £ [— /3 p ,f3 p ], e.g. eq.(B20), instead of the 
required [— /?, 0] interval where time differences appearing in propagators are evaluated (note that (3 P < (3). As 
explained in Appendix B, this departure of the standard equilibrium KMS condition in the particle number case is 
crucially different from that of a chemical potential associated to a exact charge conservation, like the electric charge 
for charged particles. In the charge case, the departure of KMS is given by the constant /iQ-dependent multiplicative 
factor in (B27l and, as explained in the Appendix, there is no obstruction to define the Matsubara representation for 
t 6 [— /?, (3\ in that case, e.g. eq.(B28l, which amounts just to a shift in the Matsubara frequencies. 

Turning again to the case analyzed in this paper, the KMS-breaking mentioned in the above paragraph may 
be a problem for instance in diagrams contributing to the partition function (closed diagrams) whenever there is 
momentum exchange (time propagation) inside the diagram, i.e, more than one interaction vertex, since in that case 
the imaginary-time variables running in the internal propagators lie in the interval [— /3, f3] while those propagators 
are only /3 p -periodic. When there is just one interaction vertex, time integration factorizes trivially and the answer 
is proportional to powers of the tadpole-like contribution Ay(r = 0,x — 0) given in Appendix |b| That will be the 
case for all the contributions to leading order 0(T 6 ) in the calculation of the ChPT partition function. The diagrams 
that contribute are given in Fig{2](see section [IV A| . However, consider for instance the diagram labeled 8b in Figj2] 
contributing to the ChPT free energy density to 0(T 8 ). Taking for simplicity constant vertices, as in the case of 
V = \<fi 4 /4l, this diagram in the ITF would be proportional to: 



G 2 (0) f d 3 p 
P J (2tt) 3 



dr' drA T (T - t',p)A t (t' - r,p) 



(20) 



with G(0) = At(0) gi ven i n eqs.(B29 l-(B31 1. Now, as commented above, we cannot just replace the Fourier rep- 
resentation for At in (B20l since it is only defined for the [— /3 P ,(3 P ) interval. This obstruction produces additional 



unnatural terms. The app e aranc e of those terms can be seen by using the mixe d re presentation for At in terms of 
G > and G < given in (B10l-(B11 1 and performing explicitly the t,t' integrals in (20). We get: 



I = 



-G 2 (0) 
G(0) 



d 
dm 2 
d 3 p 



G(0) 



(3 J (2tt)3 8££1 1 



-n p (E p )] 2 [e^ 



1] [n p (E p ) 



[e-2/3/* 



(21) 



with Ep = \p\ 2 + m 2 . The first term above gives the standard result for /i = with the replacement of the distribution 
function n — » h, as one would expect from kinetic theory arguments, while this property does not hold for the 
additional terms. The remaining contributions vanish for /i = but they do not do so in the T —t + limit where they 
diverge. This contradicts the natural physical expectation that in the T — > + limit and for /i < m, the free energy 
should reduce to the vacuum contribution. 

A related conflict arises when trying to calculate correlation functions in the ITF. The loss of KMS /^-periodicity 
implies that the dependence on external times is not only through time-differences. In particular, this means that 
correlators depend on ti. Consider for instance the tadpole-like contribution (we omit the spatial dependence for 
simplicity) : 



iti+P _ _ pT+0 

drA T (r 1 - r)A T (r - t 2 ) = / g?tA t (ti - r 2 - r)A T (r) 

iU JT 



(22) 
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~* T '-\- Q (3 T T I Q 

where T — iti — t 2 . Now, if we make J T = j Q — j Q + f„ , the change of variable t — > t + /3 in the third integral 

does not cancel the second one due to the loss of /3-periodicity of At- Therefore, the result does not depend only 
on Ti — T2 but on T, i.e., depends on t\ and t 2 independently and, as a consequence, the dependence on ti does not 
vanish. 

As we will see in sections |IV A| and |IV C| terms of the type showed above appear in the self-energy to leading order 
and in the partition function at order 0(T S ). In the latter case, our approximation reaches its validity limit, since 
particle-changing processes start playing an important role. However, precisely for that reason, at the temperatures 
where the C(T 8 ) needs to be included we may consider in practice fj, <C T,m for these contributions. Recall that 
in fact, the conflictive terms in (21 1 are O(fifT) so that we will be introducing only small corrections by neglecting 
them. The presence of those unnatural terms in the ITF may be also understood if we note that we are facing a 
nonequilibrium situation, where the ITF is not appropriate and which must be formulated using a contour including 
real times [HI [25]. We will indeed see next that one can define a suitable Real-Time Formalism (RTF) so that these 
problems are not present, at least to the order we consider here, and one can calculate properly not only thermal 
correlators but also vacuum diagrams contributing to the free energy. 



Real-time formalism 



We consider now the full contour in Fig jl] and, following the standard notation, we denote by Dij — G(ti — tj) with 
ti € Ci,tj e Cj. We then have for the C\ t % parts of the contour (we omit the spatial dependence for simplicity): 



D n (t-t') 


= G>{t- 


t')6(t-t') + G<(t 


- t')0(t' - 


t) 


£>22{t-t') 


= G<(t- 


t')e(t-t') + G>(t 


- t')0(t' - 


t) 


D 12 (t-t') 


= G<{t- 


t' + ia) = D 21 {t' - 


t) 





(23) 



where t,t' S M. and G > ,G < given in (BIO I-(B1 1 1 and so on for the remaining components. 

In order to formulate properly the RTF at // 7^ we take first, as customary, ti — > —00. This is necessary if we want 
to calculate Green functions with arbitrary real time arguments. In principle, this choice implies also that, imposing 
vanishing asymptotic conditions for the j currents and for the spectral function, which hold also in our case, the 
generating functional for V — can be factorized as [TT] : 



~Zj=« \j] = AA Z v= c ° i2 [j] Zffi M [j] , (24) 

so that one could calculate real-time correlation functions without worrying about the imaginary-leg contributions. 
However, as it was pointed out in |26l 127) . there are imaginary-time contributions that still survive in particular 
diagrams, for instance self-energy insertions, which indeed we will calculate here. Nevertheless, there is a standard 
rule for collecting all the relevant contributions but using only the propagators in C1.2, the so called |po|-P r escription 
[24), |26| 127] , This prescription amounts to use in Fourier space n(|po|) instead of the seemingly equivalent n(E p ) 



when multiplied by the on-shell (^-function, as in (B19). For instance, with this prescription one obtains that a 
simple constant tadpole-like insertion in the self energy such as the diagram showed in Fig(6^ with a constant vertex, 
amounts to a redefinition of the mass, as expected. It also guarantees that there are no ill-defined contributions, 
such as products of S distributions at the same point which in principle could appear when multiplying the RTF 
propagators. What we will show here is that for /i ^ there is also a natural prescription which works, now in terms 
of n — > n p (po), leading to the same properties at the order considered here. However, it must be pointed out that 
to higher orders, there may be additional ill-defined terms arising from a nonequilibrium distribution [25 . Our RTF 
avoids the main obstruction that we faced in the ITF, since the length (3 of the imaginary leg disappears from the 
integration limits in momentum space, whose Fourier representation is well defined now. Moreover, we also choose 
a — ► 0^ . Therefore, for Green functions with real-time arguments for which we neglect (with the above prescription) 
the 63,4 parts, we end up with a Keldysh-like contour characteristic of nonequilibrium Thermal Field Theory [16 . 
With this procedure we will see that an additional property holds: most results can be written as functionals of n, 
which encodes all the T, /1 dependence. This is also a expected property from kinetic theory arguments, at least for 
the leading order corrections in n (dilute gas regime). 

This allows then to calculate properly any real-time correlation function directly, i.e., without appealing to the 
analytic continuation from the ITF, which is cumbersome for fi ^ 0. In addition, as we will see below, one can also 



obtain information about the free energy density without using the ITF. Let us then write the propagators (23 1 in 
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momentum space for our choice of contour (note that the D\\ component corresponds to the free propagator G in 
|B19l): 



Dn(Po,p) 

£>22(P0,P) 

Dl2(P0,P) 
E>2l(P0,P) 



P0~ E p + ^ 



PI-EI + le 



2<K8(jpl-El)n{\p Q \-ix) 
27r5(p 2 ~ E 2 p )n(\p Q \ - n) 



27rS(p 2 Q ~E 2 p )[0(-p o )+n(\ Po \- f ,)} 
2ir6(p 2 -E 2 p ) [9(po)+n(\po\-n)] 



(25) 



In the above propagators, we have chosen, as discussed above, the |po| prescription ensuring that the distribution 
function does not depend explicitly on E p , as in the n = case. We will see below that this yields the same expected 
properties as for /i = 0. It can also be readily checked that our /i ^ RTF propagators above coincide with those 
given in |15j . obtained assuming a direct replacement of the distribution function by the /i ^ nonequilibrium one. In 
fact, the free propagators (25 1 can be readily recast into the general nonequilibrium Kcldy s h for m, given for instance 
in [TniUS] by taking for the nonequilibrium distribution function our n p (ko) given in (B12 1-( B13 1 3 , which satisfies the 



property (B14l, as required for general nonequilibrium derivations [25 



To provide a particularly relevant example of our previous statements, let us consider the tadpole-like correction 
to the self-energy given by the diagram in Fig[6^ with a constant vertex (the generalization to derivative vertices 
appearing in ChPT calculations will be straightforward). The external leg is fixed to be of "type 1", since we are 
calculating the two-point function with real arguments, i.e, the first-order correction to D\\. Then, if we consider 
only the C\% contributions, this diagram gives in position space: 



F(x-y) 



J'=l,2 



D lj (x-z)D JJ (0)D jl (z-y) 



(26) 



Note that correlators depend only on space and time differences in the RTF, so that the problems discussed in the 
previous section, related to the ITF version of the tadpole in eqj22J) are n ot present now. 

Now, we take into account that Dn(0) = £>22(0) = G(0) in (|B29|)-(B31 ). Then, the Fourier transform of F is: 



F( Po ,p) = iG(0) [blM.P) ~ D 12 (j) ,p)D 21 {p ,p) 
We replace in the above equation the propagators in ( 25 1 and use 



(27) 



8'(x) 



we keep the regulator in the definition of S(x) = ^ ^ x +m+ ' 



x i0+ — 2 inS (x) where as customary 

Thus, we can write: 



F( POlP ) = »G(0) 



Pi 



E 2 
v 



— 2nin 



m(\ Po \ - M )*'(pg - E 2 ) = -G(0)-^D n (p ,p) 



(28) 



Note that it is in the last step in the previous equation where i t is crucial to use the \po\ prescription chosen above 
since n(|po| — /•*) is independent of m . Therefore, the result (28 1 implies that the only modification in the Du 



propagator is m — > m — G(0), which is the expected result of mass renormalization which in addition is obtained 
from the fi = case by replacing n —* n in the (finite) thermal correction to the tadpole diagram given by the function 
gi(m,T, fi) in (B31 1. Note that to this order and with this prescription we have been able to get rid of the ill-defined 
5 2 terms. However, this prescription might not be enough when higher orders are included, since it has been shown 
in |25j for a \<fi 4 theory that additional nonequilibrium ill-defined terms arise, which should be properly regulated 
with a nonzero particle width. At the order corresponding to our previous result (28 1 we coincide with |25j . We will 
comment more about this issue in section TlV CI 

Two more important remarks are in order. The first one is that the spectral properties of the interacting theory 
are really defined from retarded Green functions, not from time-ordered ones. From the ITF, retarded correlators are 



3 The convention in 1251 is such that the D12 and D21 components are reversed with respect to ours. 
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defined directly by analytic continuation. However, we have seen that this is not a well defined procedure for fi ^ 0. 
The solution of the problem of finding retarded Green functions from the RTF time-ordered product was given in 
|28j . In that work, a set of rules (so-called circling rules) were provided in order to define a function that has the 
required causal retarded properties, namely that satisfies that one of the outgoing lines of the corresponding diagram 
has the largest time component. It was then shown in several examples that this function coincides with the analytic 
continuation of the ITF correlator. Now, it can be checked that the same properties of the free propagators used in 
[2"8] for the derivation of the circling rules hold for our propagators and therefore the same rules lead to the RTF 
retarded function at fi ^ 0. The application of those rules is trivial for the tadpole case discussed above, since there 
is only one vertex. However, they will be of use for the case of higher order contributions to the self energy which we 
will consider below, like the thermal width arising from diagram [6]}. 

The second remark has to do with the calculation of thermodynamic quantities within the RTF, i.e., the partition 
function or the free energy density. In principle, due to the factorization of the imaginary- leg commented above, the 
contribution to vacuum graphs when summing over fields of type 1 and 2 vanishes identically. However, it was shown 
in |29j that fixing one of the vertices of a vacuum diagram to be "external" of type 1 and summing over the remaining 
internal vertices with an overall p factor reproduces the free-energy result and for /i = coincides with the ITF. The 
functional arguments used in those papers are also applicable to our /i ^ case and in fact, the direct use of that 
prescription leads to the expected answers. Let us show this for the case of the 0(T 8 ) diagram 8b in Fig 2] analyzed 
in section [ill A| in the ITF. Applying the previous prescription and with constant vertices, we get that this diagram 
is now proportional to: 

zG 2 (0) ]T J D lj (x-z)D jl (z-x)=G(0)F(0) = -G 2 (0)-^G(0) (29) 



with F in (26 1. We then see that we arrive to the ITF result (21 1 but without the additional terms discussed in that 



section, since the proportionality factors between this diagram and (29 1 or (20 1 come only from combinatorics and 



are therefore identical. We will use this real-time prescription to define properly our free energy. 

IV. APPLICATIONS TO THE PION GAS 
A. Evaluation of the ChPT free energy. 

We apply our previous results to the pion gas, described by ChPT with two light quark flavors of mass m u — 
m<2 = m q |30l 131] . The lagrangian is constructed as an expansion in derivatives and pion masses, generically 0(p) 
with p < A x ~ 1 GeV, so that £ = £2 + £4 + ■ • • with £2^ = 0(p 2k ). In the range of temperatures and chemical 
potentials we are interested, both T, ^ = 0(p) formally, which corresponds to T below T c ~ 200 MeV. The ChPT 
0(p D ) power of a given diagram is given by Weinberg's power counting D = 2(N L + 1) + ^2 k 2N k (k — 1) [32] where 
Nl is the number of loops and is the number of vertices coming from £ 2 fc. In our approach, we do not perform any 
formal chiral expansion in except in higher order contributions (see our discussion below and in section IIP where 
it is reasonable to expand in fi^/T. We will follow closely the notation and conventions in [31j . where the explicit 
expressions of the £2 and £4 can be found. The lagrangian £ 2 is the nonlinear-er model, whose free parameters are 
the pion decay constant and mass to leading order / = / ff (l + 0(p 2 )) with f n ~ 93 MeV and m = m w (l + 0(p 2 )), 
TYL-K ~ 140 MeV. To fourth order, £4 contains five independent low-energy constants £i_4 and h\ which absorb the 
divergences of the one-loop diagrams with only £2 vertices. The renormalized U appear in physical processes such 
as pion scattering and therefore their values can be fitted experimentally. We will use the same central values given 
in [301 131] in order to compare more easily with the results in [3T] at fi„ = 0. Those values are l\ = —6.6, I2 = 6.2, 
I3 = 2.9, I4 = 3.5. The constant hi multiplies a contact term and appears in the vacuum free energy and quark 
condensate. We use also the estimate in [30] [31] of hi ~ 3.4. The lagrangians of higher orders will only appear 
through renormalization either of the vacuum energy or the pion mass and therefore the low-energy constants of those 
orders will not show up once the results are expressed in terms of the physical pion mass (see details below). 

The free energy density z, from which thermodynamical observables can be obtained, is defined as customary: 

z(T,n v ) = -T Jim ilogZ^T,^) (30) 

We also define the thermodynamic pressure as in |31j . i.e, subtracting its T = contribution given by the vacuum 
energy density: 
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0(T B 



0(T 8 



2 4a 4b 



(ia 



6b 6c 



° O h OO O m 

8a 8b 8c 8d 

AOOO©00 

8e 8f 8g 8h 



FIG. 2: Feynman diagrams contributing to the partition function of the pion gas up to and including 0(T 8 ). The first row 
includes diagrams up to (D(T a ) while the second and third rows are the 0(T 8 ) contributions. The dots denote interaction 
vertices coming from £2 while those vertices coming from higher order lagrangians are indicated by a square box. The notation 
is the same as in 1311. 



P(T, fa) = Z - z(T, 



zq = lim z 



(31) 



It is important to emphasize that all our results for the pressure and quantities derived from it have to be understood 
strictly as time-dependent throughout the plasma expansion, in the sense explained in section [TTJ, the time evolution 
towards a chemically equilibrated phase being driven by fi v (T). The diagrams contributing to the free energy in 
ChPT arc the closed diagrams showed in Fig(2j where we follow the same convention as [31] to name the diagrams. 
The number assigned to each diagram indicates the order in the chiral expansion and the numbers inside the boxes in 
the vertices refer to the lagrangian order, the case of £2 being indicated by a dot. Recall that for a given order of the 
lagrangian, there are vertices with arbitrary number of (even) pions due to the chiral expansion of the S'[/(2)-valued 
chiral field U — exp(i7r°r a //) where r a are the Pauli matrices and 7r a the pion field. 

The leading order S2 — —f 2 m 2 coming from the contact term (independent of the pion field) in £2, is independent 
of T and fi n and therefore contributes only to the vacuum energy density Zq. Note that, according to our discussion in 
the previous sections, we will ensure that all our contributions have a well-defined T — > + limit for \i v < m, i.e, that 
the contributions to zq to any chiral order are /^-independent. The next order corresponds to diagram 4a and 4b. 
Z4 a corresponds to the quadratic pion field contribution in £2 and is therefore nothing but the free partition function 
given in ( B34 ) multiplied by 3 accounting for the three pion degrees of freedom. The divergent contribution to z^ a is 
T and n w independent and therefore it merely renormalizes zq. 

The next order in the chiral expansion is 0(T 6 ) and the diagrams contributing are zg a bc in Fig(2j It is important 
to remark that this is the first order where pion interactions show up. Graph 6c renormalizes zg, while 6b is of the 
same form as 4a and therefore gives rise to the free partition function contribution but with the mass shifted by its 
tree level £4 renormalization (see section |IV C I which depends on ^3. As for diagram 6a, taking into account ( B32 1 



III 



its contribution is proportional to G 2 (0). As discussed in section 
both ITF and RTF and corresponds to the result in [31] replacing G(0) 



in this case the result is trivially identical in 
G(0): 



ZGa 



w 6 ^ 



(32) 



The divergent contribution in (32 1, according to ( B29[)-( B30 ), contains a contribution to z$ and another one which 
cancels, as it should, with the one in £3 so that, using (B35l, the total finite result for the pressure to 0(T 6 ) is: 



3 



3 m 2 
8/ 



2 Mm^nJY + OiT*) 



(33) 
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with the functions g\ and go given in (B31 1 and (B33l respectively and where is the physical pion mass at 



T = fi n = 0, related to the bare mass m to this order as |30j : 



m„ = m 



32tt 2 f 2 



0(m 4 



(34) 



Recall that, to this order, the difference between and m is only relevant in the go contribution in (33 1. The 
same applies to the distinction between / and / w : 



Si = / 2 



h rn 2 A . 
1 + ^ + °(™ 4 ) 



(35) 




FIG. 3: (a), (b) Diagrams contributing to leading order (tree level) to 2n — > 4n processes, (c) A one-loop contribution to 
elastic nix scattering. 

We consider now the C(T 8 ) contributions showed in Figj2j Now, there are several aspects which make the cal- 
culation qualitatively different from the 0(T e ) one. An important point is that to this order we may expect that 
our approximation of particle number conservation is less accurate, since vertices entering number-changing processes 
show up. Consider for instance the diagrams contributing to 2ir «-> 4ir processes in the thermal bath, which to leading 
order in ChPT are given by the tree level diagrams showed in Fig H)i andlJj. Now, unlike the 0(T e ) case, one can 
draw vacuum diagrams from these processes by identifying external lines. For instance, joining lines in pairs in the 
graph in Fig(3^ as 1-2, 3-4, 5-6 and equivalent combinations leads to diagram 8a in Fig(2] Similarly, joining 1-2, 3-6, 
4-5 in Fig|3j) produces diagram 8c. This is not a one-to-one correspondence. For instance, joining 1-3 and 2-4 lines 
in the elastic one-loop diagram in Fig(3j: yields also diagram 8c. Diagram 8b can also be obtained from a elastic 
process (Figj3b joining 1-2, 3-4) or from an inelastic one (Figj3j) joining 1-3, 2-6, 4-5). The crucial point is that none 
of the 0(T 6 ) vacuum closed diagrams in Fig(2]can be obtained from the lowest order inelastic diagrams in Fig|3]i|3|3. 
This distinctive feature can be interpreted as a way to identify the validity range of our approximation. However, 
we should bear in mind that these 0(T 8 ) corrections are meant to be relevant only very near T c [31] and therefore 
in the region where chemical equilibrium is nearly restored and n„ — > 0, not surprisingly due to the presence of 
the particle-changing processes just discussed [8]. Precisely for this reason, the /^-dependence of these diagrams is 
suppressed in powers of n^/T and n^/m^. Therefore, numerically our approach will be still justified to this order. In 
addition, as we have explained in section |5l] taking fi^/T small justifies in practice to get rid of unnatural terms in 
the ITF formulation. 

With the above considerations in mind, we proceed to evaluate the 0(T 8 ) diagrams in Figj2] Graph 8h renormalizes 
So and graphs 8f,g renormalize the pion mass to 0(m 6 ). Graph 8a is proportional to a third power of the propagator 
at the origin, with the same coefficient as in |31j : 



Z8a 



48F G (0) 



(36) 



which contains divergent contributions, according to (B29l 



The graph 8b has been analyzed in section III The relevant integral contributing to this graph is ( 20 1 in the ITF 



and (29) in the RTF with the prescription discussed in that section. The difference between both formulations is of 



0(fi^jT) and therefore expected to be numerically small, for the reasons just discussed. The rest of the contributions 



13 



to this graph are proportional to G 3 (0) and the proportionality constants are the same as in [3T]. Thus, adopting the 
RTF prescription, we get: 



zsb 



16/ 4 



G^(0) 8 + 3m z 



d 
dm 2 



G(0) 



(37) 



whose divergent contribution can be also separated using (B29l 



Graphs 8d,e have the same form as graph 6a in (32 1, but due to the form of the £4 lagrangian and following also 



our previous RTF prescriptions, we arrive to the same structure as in [3Tj : 



e = --^ j(2; 1 + 4Z 2 ) 



G 



fJLU 



G(0) 



3^6 



(3Zi +l 2 + k) m 4 G(0) - 



dm 2 



G(0) 



(38) 



where G^„ = d^d v G(§), which has the same properties as in |31j , namely, its divergent contribution is the same (T 
and [in independent) while its finite part can be written in the same way in terms of go and g\ . 
The remaining graph is 8c. Following again the RTF prescription, this contribution is: 



where: 



%8c — 



48/ 4 



3m 4 Ji - 72 J 2 + 16m z G(0) 



(39) 



Ji = i / d x 



= I d 3 x I dt 



G>(t,x) - G<(t,x) 



J% = i I d A x 



(d„D n (x)d»D n (x)y - (d„D 12 (x)d»D 12 (x) 



(40) 



Written in the above form, it is not difficult to show that for fi^ = 0, when the propagators are /3-periodic, i.e. 
G < (t + if3) = G > (t), one has for instance J\ = J d 3 x drA^r, x) and similarly for J 2 [3T]. As we have seen, for 
/ijr 7^ the perio dici ty condition does not hold. However, for this diagram, instead of working directly with the 
RTF expressions (40 1, we will make use of the fact that G < (t + if3) = G > (i) + 0(13^^) and neglect the non-periodic 

terms, so that we end up with Ji ~ J d 3 x drA|,(r, x) and J 2 — J d 3 x dr ^9 m At(t, x )9 m At(t, x)^j . This 

approximation simplifies considerably the renormalization of this graph, since now we can follow th e sam e steps as in 
[3T] . First we separate A(t,x) = G > (-ir,x) = A(t,x) t =^=° + 5A(t,x) using the representation ( |B10[ )-( |Bll| . The 
divergent contributions in the integrals (40 1 are then contained in the (SAA ) 2 , 5A(A ) 3 and (A ) 4 terms and can 
be renormalized with the same counterterms as in |31j replacing the go t i by go,i- The finite part of the Ji ]2 integrals 
can be evaluated numerically. A crucial point is that this approximation is consistent, as far as renormalization is 
concerned, with our previous evaluation of the .Ssafcde diagrams since the divergent parts of the terms proportional to 
5oiffo5i an d g\ arising from the Ji. 2 integrals cancel exactly with those coming from the other four diagrams, while 
the terms proportional to g\ add together to renormalize the physical pion mass according to the definition 



m n = - lim T log P(T, u n = 0) (41) 

T-»0+ 

In addition, as it happens for [i v = 0, this ensures that neither the tree level constants from Cq nor the T, [i^ 
independent renormalization constants needed to render Ji j2 finite appear in the final expression for the free energy 
once it is expressed in terms of . We remark that with our representation, not only the renormalization procedure 
is consistent, but the final answer for the full 0(T S ) contribution amounts to replace n(E p ) — > n(E p — fi^) in all the 
spatial momentum integrals, without dealing with unnatural terms, like those discussed in section [ITT] 

After the previous detailed evaluation, we arrive finally to a finite expression for the free energy, suitable for 
numerical evaluation, with the approximations discussed above implying that the 0(T 8 ) corrections are reliable only 
for small /i^. From this expression we proceed to present our results for the /i^ dependence of several relevant 
observables. 
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FIG. 4: Results for the pressure, quark condensate and entropy over density ratio at different chemical potentials and to 
different orders in the ChPT interactions. 



B. Results for thermodynamical observables 



From the energy density we obtain the quark condensate (the order parameter of the chiral transition), the entropy 
density and the pion number density in the standard way: 



(qq) (T, fa) 

s(T, fa) 
h(T, fa) 

.2 f2\ 



(99) (0,0) 

dP(T, fa) 

dT 
dP{T, fa) 

dfa 



1 



c dP(T,fa) 
f 2 dm 2 



(42) 

(43) 
(44) 



where c = 1 - m 2 (4/n + l 3 - l)/(327r 2 / 2 ) + 0(m 4 ). 

We plot our results in FigH The first feature we observe is that the C(T 6 ) and the ideal gas curves are very 
close to one another for all the range of temperatures and chemical potentials showed. Sizable differences due to the 
interactions only show up numerically when including the 0(T 8 ). This is a also a feature of the fa — calculation 
[3"T] . For instance, in the chiral limit (m ff = 0) and for fa = 0, the 0(T 6 ) in ( |33| vanishes identically, while the 
0(T 8 ) survives, producing conformally anomalous contributions to the pressure [33]. In FigjZJwe also compare our 
results with the virial gas approach |14j , where the pressure can be written at low pion density in terms of the pion 
scattering phase shifts. In the curves snowed in Figji] the phase shifts have been calculated perturbatively to 0(p A ) 
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T(MeV) 

FIG. 5: Dependence of Hn(T) in the isentropic approximation, with the fixed value s/n — 4. 



in ChPT and using the same set of low-energy constants that for our perturbative results with the approach of the 
present paper. We see that our 0(T S ) results with p n ^ lie reasonably close to the virial result, at least for not 
very high p,^. This is a good consistency check of our present approach. 

Another general feature that we observe in the curves is that the effect of the pion chemical potential is always to 
increase thermal effects. Effectively, it acts as similarly as a reduction of the effective pion mass (this is more accurate 
for T <C m w where the typical momenta in the distribution functions are p — 0(WT/m w )) and therefore for fixed 
T, the results for increasing p v go qualitatively in the same direction as for increasing T with p v = 0. For instance, 
we see that the pressure increases for increasing p,^ and approaches faster with T the asymptotic limit P ~ 7r 2 T 4 /30 
|31j expected in the chiral limit (T 3> m n ,p w ). The effect of interactions is also to increase the pressure, producing 
additive contributions in the ChPT expansion. 

The curves for the quark condensate show that the chiral restoration temperature goes down for \i v ^ 0. This 
is also a consequence of the above discussed qualitative behavior, since the system for p v ^ is closer to chiral 
restoration. With the numerical values we get, we see that if chemical freeze-out takes place for temperatures below 
the chiral phase transition, then we do not expect to see any change in the value of T c . On the contrary if T c ^ em > T c 
(which is less likely with the available experimental information) we would expect a reduction in T c compared with 
the estimates taking p„ = 0. 

It becomes clear from our discussion in section [TT] that incorporating additional physical requirements allowing to 
describe ^(T) is crucial in our approach, in order to be consistent with the chemical nonequilibrium evolution. In this 
sense, a very interesting observable is the ratio of entropy density to pion density, also plotted in Fig(4] It has been 
pointed out [21 13 [8] that on general grounds one expects this ratio to remain almost constant during the expansion. 
This is the isentropic expansion approximation, which is exact in the high T limit T 3> m ff , p v for the ideal gas. We 
remark that we are restricting here to the gas of pions. If heavier degrees of freedom arc included, such as the p, one 
has to account for the total number of pions = n w + 2n p + . . . which includes those "stored" in the p if the channel 
p — ► 7T7T is considered as the only source of pion number changing, and similarly with other resonances (see details 
in [HIH]). The idea is then that by fixing s/n to a given value at the chemical freeze out temperature T c ^ eTO where 
/Jtt = 0, going down in the temperature scale one can keep s/n fixed by increasing p,„, as it can be seen in Figg 
This provides the isentropic dependence u^(T), which is given in [2, 8 in the ideal gas approximation. We plot in 
Fig j5] the isentropic curves ^(T) with a reference value s/n = 4, for which T c h em — 190 MeV for the ideal gas. The 
obtained curves follow a roughly linear behaviour, as expected phcnomcnologically |3J. The most significant effect we 
observe is the reduction of T c h em when 0(T S ) or virial interactions are included. This is a very natural effect since, 
as we have discussed in previous sections, that order in the interaction is the one where particle-changing processes 
begin to be relevant and drive the system back to chemical equilibration. The virial curve lies reasonably close to our 
perturbative 0(T S ) since the two approaches differ significantly only for rather high values of p n and T, which are 
not reached along the curve /^(T). In fact, in the isentropic evolution our C(T 8 ) approach is better justified since 
Ptt(T) <C T, m w . We also remark that the same effect of faster equilibration is seen when comparing the curves of 
the ideal pion gas with that of the ideal pion+resonances gas, as done in [2j. One can check that the curves in that 
paper for s/n as a function of T for different p,^ are systematically lower when including resonances, as in our case in 
Figfflwhen including the 0(T 8 ) or in the virial case and therefore the free pion and resonance gas equilibrates faster, 
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which is the feature that we are able to reproduce here including higher order pion interactions. 



C. Self energy: pion thermal mass and width. 



(a) 




FIG. 6: Diagrams contributing to leading order to the real (a) and imaginary (b) parts of the self-energy. 



Within the real-time formalism developed in section [Iff B[ we can calculate the pion self-energy for fi n ^ 0, whose 
leading order corrections to its real and imaginary parts are given by the diagrams in Figjf^i andjBJj, respectively, all 
vertices coming from the C2 lagrangian. 

It is important to remark that when nonequilibrium distributions are considered, as it is our case here, it has been 
pointed out that additional <5 2 -like or pinching-polc ill-dchncd contributions arise 25j, which should be regularized 
keeping a nonzero particle width. We will discuss the role of those contributions in the last part of this section. 

Consider first diagram [(3^. It includes a contribution with a constant vertex proportional to m 2 G(Q)dD\i(p) / 'dm 2 
which directly renormalizes the pion mass, following the prescriptions explained in section [lII B| and derivative vertices, 
which contribute either proportional to □G(0) = — m 2 G(0) (mass renormalization) or as G(0)p 2 dDn(p) / dm 2 — 

G(0) (Du(p) + m 2 dD\i(p) / dm 2 ^ (mass and wave function renormalization). One can then follow similar steps as 

in the standard derivation of the thermal corrections to the pion self-energy to this order [mi EH], the wave function 
renormalization being directly related to the thermal / w through the usual definition in terms of the residue of the axial- 
axial current correlator. The ultraviolet divergences arising in the calculation are absorbed by the renormalization of 
the low-energy constants I3 and I4. We finally obtain: 



fl{T,^) = f 2 -2~ 9l (m,T,^) + 0(m 4 ) 



(45) 
(46) 



with niir and /„■ the T = fi n = physical values given in (34 1 and ( 35 1 in terms of m and / to this order. 

Taking into account now the corrections to the quark condensate to the same chiral order, i.e., 0(T 6 ), which is 
given from ( 42 ) and ( 33 1 using ( B35 ) : 



(qq)(T,^) = (mXO.O) 



1 - ^p9x{m,T,^) 



0(T 8 ) 



(47) 



we obtain that the Gell-Mann-Oakes-Renner (GOR) relation [35] holds also for ^ to this order (one-loop ChPT): 



f 2 (T,^)ml(T,^) _ /2(0, 0)774(0,0) 



(qq) (0,0) 



(48) 



The GOR relation in terms of thermal quantities at T 0, fi„ = had been verified to one-loop in [35]. To this 
order, the thermal mass varies little, also at fi n ^ (see below) so that the evolution of f n follows that of the quark 
condensate and both behave as order parameters. However, beyond one-loop, the GOR does not hold for thermal 
quantities [36] . 

Another important observation is that the shift ( 45 1 in the mass to this order can be written, as in the = case 



[10] in terms of the elastic pion-pion forward scattering amplitude, from (B31l: 
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ml(T,^)-ml = -J ^^ -^J- y — Ite [T4( S = (£ p + m OT ) 2 -|/f)] (49) 
where E 2 — m 2 + \p\ 2 and T/ ff (s) is the isospin averaged forward scattering amplitude: 

1 2 2 °° 

= r ff7r ( S ,0 ! u) = -^(2/+l)T 7 ( S ,0,n) = — ^^(27 + l)(2J + l)M«) 

7=0 7=0 J=0 

^[too(s) + 9t 11 ( S )+5t 20 ( S )] = -^ + 0( S 2 ,m 4 ) (50) 

where the last expression is the lowest order 0(p 2 ) (tree level diagrams with £ 2 vertices), Tj(s, t, u) are the projections 
of the scattering amplitude with definite isospin /, s,t,u are the Mandelstam variables satisfying s + t + u = 4m 2 , 
and tu are the partial waves, defined in the center of mass frame with definite isospin I and angular momentum J. 
We follow the conventions in |30j . In the previous expression, we have included only the partial waves with lowest 
angular momentum J < 1. Those with J > 1 are negligible for y/s below inelastic thresholds, such as the KK one, 
and for the temperatures involved here [10] . 



The result in ( 49 1 is the generalization to fi^ ^ of the formula relating the shift in the self-energy with the density 
of states and the scattering amplitude to lowest order in the density (dilute gas regime) [10] . These type of relations 
were first derived by Luscher |37j studying finite- volume corrections. A very interesting point is that it admits a 
natural extension |10j by considering (in the dilute gas regime) not only the perturbative tree-level 0{p 2 ) amplitude, 
but also higher orders, including unitarized amplitudes. In the latter case, unitarized partial waves tYj(s) for irir 
scattering can be constructed to satisfy exactly the unitarity relation: 



lratUs) = ^l~^\tU (51) 

matching at the same time the perturbative ChPT expansion and providing expressions that can be analytically 
continued to the complex s plane. All these features are satisfied by the Inverse Amplitude Method scattering 
amplitudes |38j which reproduce scattering data up to \fs ~ 1 GeV and all the low-lying resonances, which in the 
pure pion case considered here reduce to the p(770) and the /o(600) or a. Recall that the ChPT amplitudes satisfy 



the unitarity relation (511 only perturbatively order by order, violating the unitarity bounds for higher energies and 
thus not being able to reproduce resonances. 

In FigjTJwe show our results for the thermal mass, considering 0{p 2 ), 0(p 4 ) and IAM unitarized amplitudes in (49 1 



We have used the same set of low-energy constants as in our previous calculations, i.e, the U given at the beginning of 
section [IV A[ For the case of the unitarized amplitudes, this set is adequate to compare with the perturbative ChPT 
expressions, although it does not give the best results for the mass and width of the resonances generated with the 
IAM. We have checked that our results do not change qualitatively by changing for instance to the set given in |39j . 
which gives better physical values for the p, /o(600) mass and width. 

Our results show that the leading order, the ChPT 0(p 2 ) given by the tadpole diagram in FigRm, produces a 
thermal mass slightly increasing with temperature and chemical potential. However, including the 0{jr) or unitarized 
corrections to the amplitude, the mass tends to decrease significantly with T and /i^. Our results at \x v = agree 
with [10 . The difference between the 0(p 4 ) and the unitarized curve is not very relevant here. Our 0(p 2 ) curve 
agrees reasonably with a linear-sigma model calculation [40] which agrees with ChPT to this order at \i v = and 
where the chemical potential is introduced in analogy with the charged scalar field case. 

These results suggest an interesting scenario: the pion system could undergo Bose-Einstein (BE) condensation 
driven by the dropping of the thermal mass by interactions. Recall that we are dealing with BE condensation of 
both neutral and charged pions, since we are considering an electrically neutral system with finite pion density. The 
physical situation is then different from the charged pion or kaon condensation taking place in nuclear matter |41j or 
isospin chemical potential [5D] scenarios, although the dropping of the effective mass takes place also in the former. 
BE condensation for pion number and its possible phenomenological consequences in heavy ion collisions has been 
extensively studied in the literature [42 . Among the observable consequences are the anomalous enhancement of the 
low-pT pion spectrum and of number fluctuations in high multiplicity events. 

In our grand-canonical interacting framework we can describe the corrections to BE condensation due to pion 
interactions. In the standard free case, the BE limit is reached when fj,„ — > m, from below (by definition the system 
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FIG. 7: Results for the thermal mass dependence on temperature and pion chemical potential, considering different orders in 
the scattering amplitude in (491. 



is below the condensed phase). Those values for the pion chemical potential seem too high compared with those 
measured in heavy ion collisions at thermal freeze out T t her ~ 100 MeV [H In other words, the required 
pion densities for BE condensation might not be reached. However, if the effective particle mass m ff (T, drops by 
interactions among the medium constituents, the value fi n — m 7r (T, /i^) would be reduced. We show that line in Figj7] 
In Fig§(left) the resulting fi BE (T) = m^{T,^ E (T)) curve is represented and compared with the isentropic curves 
corresponding to different values of s/n. We see that the BE curve thus defined lies not very far from the isentropic 
approach and the expected phenomcnological values. Those curves correspond to the 0(p A ) amplitudes, both for the 
thermal mass and for s and n (in the virial approach). In Fig[8] (right) we show also the density-temperature curves 
corresponding to ^ — > m~ for different orders in the interaction. The 0(T e ) allows for lower density values, but the 
virial contribution points in the opposite direction. We also show the curve corresponding to the BE limit by lowering 
of the mass, as we have just explained, for the same virial approach, which produces a considerable lowering of the 
required densities. In any case, the corrections due to interactions are small near thermal freeze-out. We also remark 
that some of our previous results, including those regarding BE condensation by mass reduction, rely on the validity 



of the dilute gas regime, for instance when using (49), but corrections might be important for temperatures close to 
the chiral transition or chemical potentials close to m w . 

Finally, we turn to the calculation of the leading order imaginary part of the pion self-energy, given in 
ChPT by the diagram in Fig|6|}. This is the leading order contribution to the thermal collisional width T p — 
— Im£fl(E'p, \p\)/(2E p ) <C E p with E p = \/\p\ 2 + m% and the retarded self-energy, which defines the dispersion 
relation p 2 = ml + E fl (p ,\p\;T). 

As we have commented in section |llI B| we will evaluate the retarded correlator in the real-time formalism, following 
the circling rules in |28j . which also apply to the fi n ^ case. Applying those rules to diagram |6|d we have: 



1 



2 



1mZ R (p o ,\0\) = -- H>( Po ,\p\)~H<(p 0l \p\) =--[E 2 i(po,|£|)-Eia(PD,|fl)] (52) 



2 



where H > ^ are obtained by using for the three internal lines in the diagram the G > ' < ' RTF propagators. With the 



usual RT self-energy definition [24] and our convention for the -D12, -D21 propagators given in section IIIB we have 
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FIG. 8: Bose-Einstein condensation lines. Left: the curve fj,^ E = m-niT, /j,^) with the thermal mass from the 0(p 4 ) amplitudes, 
compared to the isentropic expansion curves for the virial case to the same order and for different s/n values. Right: pion 
density versus temperature in the BE limit ^ — > m~ for different orders in the interaction, compared to the ideal gas and to 
the virial case with thermal mass. 

H > = iS 2 i and H < = i£ 12 for diagram 6b, since, once a particular choice of ij indices (i,j = 1,2) has been made 
for the two vertices in that diagram, the three internal lines carry the same ij combination. 



Now, according to our discussion in section III B and in Appendix [B 



G>(k) = 2n5(pl - E 2 p ) [0(p„) + n(|po| - /**)] = e^ k °G<(k) = e 0[*o-/M»gn(*o)]6<( fc ) 

(53) 

so that we get for the thermal width: 

r p (T, Mff ) = ^jf[^A(k 1 ,k 2 ,k 3 , P )G>(-k 1 )G>(k 2 )G>(k 3 ) 

x (l - e -^-^f(^,k 2 M))\ (2ir) 4 S(E p + jfe? - k% - k°)5 (3 \p + jfei - - fc 3 ) (54) 
where k\ 2.3 label the three internal lines, A is the squared vertex function coming from the £2 lagrangian and: 



/(*!, k 2 , ka) = sgn(k° 2 ) + sgn(fc 3 °) - sgn(fc°) 



(55) 



Recall that in the p^ = case, the /-term is absent so that one ends up with a prefactor e~@ Ep — 1 = 1/ (l+n(E p )) in 
the thermal width. The natural expectation from replacing just the distribution function n — * h p for ^ would be 
then E p — *■ Ep — fi^ in that factor, as well as the modifications of the internal distribution functions n(Ei) — > n(Ei — fj,^) 
where Ei is short for E^. This is indeed the result found in [§] derived from kinetic theory. In our case, it is not 
obvious that the answer is the same, since the function / above is not equal to one for the eight possible combinations 
of signs of the three internal We denote them by S1S2S3, with s$ = sgn(fc°). Now, we take into account that the S 
functions in each of the internal lines put them on-shell, i.e., fc° = ±Ei and global energy-momentum conservation in 

the diagram imposed by the 8- functions in (54 1. Thus, the combination H giving / = — 3 is excluded by energy 

Ei > > — E% — E3. On the other hand, from three-momentum conservation and the on-shell 

E 2 p + El -El -El and 

-C < B < C with C = 2(pki + fefe) and where hi, p are short for \ki\ and 
— (/ = — 1) is also excluded, since for that combination E p + Ei = E 2 — E3 
The same reason excludes H \- (/ = — 1). Combinations — h + (/ = 3) 



conservation E„ 



u p -r ^1 ^ yj j> —rj2 — -^3 
conditions we have that for any combination it should hold A = B where we denote A 



B — 2(fc 2 • &3 — p ■ k{) and, in addition 
\p\ respectively. Therefore, the case + + 
so that A = -2(E p Ei + E 2 E 3 ) < -C. 

and (/ = -1) give A = 2(E p E 1 ■+ 

remaining arc + + +, — I — and h 

as for /Ltjr = 0, given in [10] . 



E2E3) > C and are thus excluded as well. Therefore, the only combinations 
, the three of them giving / = 1 and the same contribution from the vertices 
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It is not difficult to repeat the above analysis, now with the external energy po — —E p . In that case, every 
combination of relative signs between the Ei is obtained from the previous case by flipping the three Si, the A, B, C 
functions being independent of the sign of Pq. Thus, applying the same arguments, the only surviving combinations 

are now , H h and + H — , the three of them giving / — — 1. Therefore, what we have proven in terms of the 

12 and 21 components of the self-energy for this diagram is: 

Eg(po = ±E p ,p) = e ^°-^ s S n fr°))£f*(p = ±E p ,p) = eM^ po = ±E p ,p) (56) 

which is the usual equilibrium relation with [3 — > f3 p . This relation will be of use in the discussion at the end of this 
section about possible higher order corrections related to pinching poles. 



In conclusion, the result (54 1 we find with our diagrammatic method is the same as in kinetic theory [9], which, 

after relabeling k\ <-> — k 3 in — I — and k\ «-> — k 2 in h and performing the three integrals in kf using the on-shell 

(^-functions, can be written as: 

I i r JL ^3 . 

F P (T,^) = -jj, r ]J w ^r-n(E 1 -^)[l + n(E 2 ~^)}[l + n(E 3 -^)} 

8E p 1 + n(E p - (in) J fJi (2iry2Ei 

x \T^(s,t)\ 2 (27r) 4 S(E p + E 1 -E 2 -E 3 )S^(p + k 1 -k 2 -k 3 ) (57) 
where is the isospin averaged elastic pion scattering amplitude with s = (E p + Ei) 2 — \p+ k\\ 2 , t = (E p — E 2 ) 2 — 

Taking now the dilute gas regime in the previous expression, which amounts to neglect all the Bosc- Einstein functions 
n <SC 1 except n(E\ — fj,„), gives rise to the extension of Luscher's formula in terms of the forward scattering amplitude, 



as in (49 1 but now for the imaginary part of the self-energy through the pion thermal width (which vanishes at T = 0): 



rDG / T \ 1 f d3 k (T? , Vs(s ~ 4m 2 ) If d 3 k f 

^ = 2E p J (2nf n{Ek - ^ 2E k ^ (s) = 2E p J J^f2E k n{Ek ^ ImT ™^ 



(58) 



where we have relabeled k\ —>■ k and ov^ is the total irir cross section 



<r*n(s) = ^y;(2j+i)(2j+i)|t /J ( S ) 
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, , 1 A 2 , Imr/ ff ( a ) (59) 
y s(s — 4m z ) 



where the last line is the optical theorem, following from exact unitarity (51 1. 

We remark that our final results both for the real and imaginary parts of the self-energy to this order correspond to 
the replacement n —* h evaluated at positive energies. This is not only natural from the kinetic theory viewpoint but 
it is also formally obtained by performing such replacement in the analytically continued fi^ = ITF self-energies. 

The thermal width is of phenomenological relevance, since it enters directly in the calculation of transport coefficients 
in the meson gas |33L 143] . It is then important to estimate pion chemical potential effects in the width during the 
phase of chemical nonequilibrium where particle number is approximately conserved and transport phenomena can 
be described relying on the dominance of clastic collisions, which is also consistent with the dilute gas regime. On the 
other hand, in this regime the mean collision time defined for ultrarelativistic particles as r — l/(2f) [3, 8, 9 with 
the averaged width: 

fV- \ J d3 P T p (T^ n )n(E p - 

r(T,M^) = - r ,„ r (60) 

J d 6 p n(E p - 

provides direct information about thermal relaxation. We represent r in the dilute approach in Fig(9] using the 



scattering amplitude in (58 1 to different orders, including the unitarized case. We use the same set of U constants as 
in the rest of the paper. 

We see in the figure that the effects of reproducing correctly the energy behavior of the scattering amplitude is 
important for evaluating the collision time. In particular, the unitarized curve shows important differences with the 
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FIG. 9: Mean collision time in the elastic and dilute limit, considering different orders for the pion scattering amplitude and 
different values for the pion chemical potential. 



perturbative ones in the temperature range showed. This was also noticed in [10] at fj,„ = and the importance 
of including unitarized corrections to the width for transport coefficients in the meson gas has been highlighted in 
[33U43] for instance regarding violations of AdS/CFT bounds for the shear viscosity over entropy ratio or correlations 
between the bulk viscosity and the conformal anomaly. 

Another clear effect that we observe is the reduction of the mean time with the pion chemical potential, also observed 
in [9 with 0(p 2 ) amplitudes. Physically, in the temperature regime where r is much smaller than the typical plasma 
lifetime (~ 10 fm/c) which at the same time is small compared to the inelastic collision time driving the system to 
chemical equilibrium, the system remains in thermal but not chemical equilibrium. From the estimates of the inelastic 
collision rates given for instance in [5j and the results in Figj9] this would happen at fi v = typically in the range 120 
MeV< T < 180 MeV. However, precisely in that regime, and as we have explained in this paper, ^ 0, its typical 
values being given by the isentropic curve in Fig|5] which means that the range of thermal equilibrium enlarges 
from below, from the commented reduction of r with fi w . In fact, estimating the thermal freeze out temperature 
T t ) ler as that where this approximation ceases to be valid, i.e., where r ~ 10 fm/c (this type of dynamical condition 
has been used also in [T] to determine the freeze-out conditions) gives a shift in the thermal freeze-out temperature 
ATther — ~20 MeV with respect to the fi^ — case, following approximately the isentropic values in FigjEJ In 
particular, using the unitarized results in Fig|9]we obtain in this way T t her — 95 MeV, close to experimental values. 

Finally, let us comment on the pathological nonequilibrium terms found in |25j in a g 2 </> 4 context. Those terms are 
of the type of 5 2 functions at the same point, or pinching singularities and therefore have to be regularized by keeping 
a nonzero particle width in the propagators, i.e., T p ^ in our case. The first nonvanishing term of this kind in the 
g 2 </> 4 theory is the three-loop diagram given in Figure 2b of [25]. Note that formally this is an 0(g e ) correction, while 
the diagrams we have considered here in FigjoJ for which there arc no such pathologies, would be 0(g 2 ) and C(<? 4 ) 
respectively in that counting. Nevertheless, the argument in [25] is that those contributions are proportional to the 
inverse width 1/T with T = 0(g 4 ) and therefore could become of the same order as the leading ones. The form of 
such leading pinching-pole term [25 in our case is: 



dPA R (p)A A ( P ) {[(1 + n p ( Po )] £?!j - Yf&nM) (61) 

where the self-energy components of diagram 6b that we have analyzed above enter directly to this order and A^a 
denote the retarded/advanced propagators. Let us isolate the leading T p — > + behaviour of the previous expression. 
The product ArAa is the characteristic pinching-pole contribution appearing typically in diagrammatic calculations 
of transport coefficients [33J 113] and in the T p — > + behaves as: 



A R (p)A A (p) 



5(P 2 -E 2 p ) 



(62) 



which puts on-shcll (po = ±E p ) the integrand of (61 1. Now, we recall the relation (56) we have derived for diagram 6b, 

together with the properties of the "modified" distribution f unc tion, in particular 1 +n p (p ) = s l3p n p (po). Altogether, 
this means that the leading order pathological contribution (61 1 vanishes in our case in the T p — > + limit. There may 



be higher order corrections of this kind, but in accordance with the power counting in 25j, those would be subleading 
with respect to the contributions in Fig(6]that we have analyzed in this section. 

The previous analysis showing the absence of pathological terms is only valid to lowest order in those terms. Further 
conclusions can only be reached with a complete ChPT analysis of higher order pinching-pole contributions, extending 
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that in j^S], which is beyond the scope of this work. Nevertheless, it is worth pointing out that the low-T ChPT 
counting does not involve any coupling constant, which implies important differences with respect to the g 2 4> 4 one 
in the perturbative behaviour of pinching diagrams at low T |43j . In addition, the analysis of [25] shows that the 
nonequilibrium pathological terms are always proportional to 8n, the deviation from the Bose-Einstein distribution 
function. Thus, in our case, we can use an argument similar to the one we invoked in section IV A when dealing 
with the C(T 8 ) terms in the partition function. Namely, that in the ChPT counting those diagrams are expected 
to be important for temperatures for which /3/i(T) ~ Sn <C 1, giving a further suppression. We finally remark that 
in the fx n ^ scenario, the presence of non vanishing pinching-pole terms have also their origin in the presence of 
particle-changing processes (see also our discussion in section IV A I. In fact , self-energy combinations of the form ( 6 1 1 
enter directly in the Boltzmann equation describing the rate of particle number change |25j . The fact that the leading 
correction of that type (61 1 vanishes in our case seems to be related to the fact that the leading self-energy corrections 
arising from the diagrams in Fig(6] can always be expressed in terms only of the clastic nir scattering amplitude, as 
we have discussed extensively throughout this section. 



V. CONCLUSIONS 



In this work we have developed a path-integral diagrammatic formalism in order to deal with chemical nonequi- 
librium effects in interacting scalar field theories, in the regime where particle number is approximately conserved. 
Within the theoretical framework of holomorphic path integrals and thermal field theory, we have derived the rele- 
vant Feynman rules for nonzero particle number chemical potential /i, whose validity is restricted to the temperature 
regimes where one can neglect particle-changing processes. This derivation in the interacting case is original of this 
paper, to the best of our knowledge. 

We have addressed some subtleties related to the choice of contour in complex times, leading to the extension of 
real and imaginary time formalisms at [i ^ 0. We have shown that the consistent formulation is the real-time one, in 
agreement with other nonequilibrium formulations. The imaginary-time formalism can lead to spurious contributions, 
related to the loss of periodicity or global KMS conditions. These problems are not present in the real-time formalism, 
once a proper energy representation for the propagators is chosen, in accordance with the standard /i = choice. In 
addition, following previous studies in the literature at \x — 0, we have been able to construct the combinations of 
real-time diagrams leading to retarded correlators and to closed diagrams contributing to the free energy. 

We have applied this formalism to the case of a pion gas, relevant for Relativistic Heavy Ion Collisions between 
thermal and chemical freeze-out with nonzero pion number chemical potential /z ff (T). Our description is consistent if 
the T-dependence of \i v encodes the time evolution of the plasma between those phases. The relevant diagrammatic 
scheme for temperatures below chiral restoration is Chiral Perturbation Theory. In this framework, we have calculated 
the leading corrections to the ideal gas coming from chiral interactions. To leading order C(T 6 ) the corrections to 
the pressure can be expressed in terms of tadpole diagrams and are numerically rather small up to T c . To next 
to leading order 0(T 8 ), closed diagrams contributing to the free energy can be obtained from particle-changing 
processes, which signals the onset of the number conservation approximation break-up. Nevertheless, since //^ is 
small for temperatures where those ChPT corrections become important, they can be reliably calculated and produce 
sizable deviations from the free gas. The results to that order agree reasonably well with a virial expansion analysis. 
Our results for thermodynamical observables show that both chiral interactions and [i^ tend to increase the pressure. 
The chiral restoration critical temperature decreases with increasing /i^, which would be of relevance only if chiral 
restoration takes place for lower temperatures than chemical freeze-out. We have also calculated the isentropic ^{T) 
curve for different orders in the interactions. The corrections to the ideal gas show a significant reduction of the 
chemical freeze-out temperature, which is the expected effect of interactions, since they increase the probability of 
producing inelastic processes. The same effect had been observed previously in a free gas of pions and resonances. 

Our approach allows to derive thermal corrections to the pion self-energy at fi^ ^ 0, from the leading-order 
ChPT diagrams, both for the real and imaginary parts of the retarded correlator. The imaginary part comes from 
a two-loop diagram, for which the use of RTF rules for the construction of the retarded function is crucial. After a 
detailed evaluation, our diagrammatic result is shown to coincide with the expected expressions from kinetic theory 
arguments. We have also discussed the role of higher-order pinching pole contributions to self-energies, providing 
different arguments which support that those corrections are subleading in our approach. The real part of the self- 
energy gives the thermal mass, which together with the condensate and the pion decay constant to the same order, 
satisfy the /i,^0 extension of the Gell-Mann-Oakes-Renner relation. In addition, both the real and imaginary parts 
satisfy a Luscher-like relation in terms of the forward pion scattering amplitude. This relation allows to calculate in the 
dilute regime the self-energy corrections for higher orders in the ChPT amplitudes, including unitarized expressions 
which have the physically expected energy behavior and reproduce the lightest resonance states. The results for 
the thermal mass show a clear decreasing both with T and for 0(p 4 ) and unitarized amplitudes. This suggests 
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the interesting possibility of reaching Bose-Einstein condensation when the effective thermal mass approaches the 
chemical potential. This mechanism would require lower pion densities to reach BE condensation. We have discussed 
this possibility, which is a purely interacting effect, within the isentropic values and comparing the pion densities 
with those in the standard approach of considering the ideal gas BE limit /i — > m~ with the vacuum mass. Our 
mass-dropping BE curve is not far, but still above the isentropic ones for reasonable values of chemical freeze-out. 
Finally, using also the scattering amplitudes, we have evaluated the corrections to the mean collision time at /i^ ^ 0. 
The mean time decreases with T and fi^ for all orders in the interaction, which implies a sizable reduction, compared 
to the fx v = case, of the thermal freeze-out temperature, estimated as that where r equals the typical plasma 
lifetime. 

Summarizing, the diagrammatic field-theory scheme developed in the present work provides, in our opinion, useful 
results regarding the chemically non equilibrated phase of the meson gas resulting from a Relativistic Heavy Ion 
Collision. In future works we plan to generalize the analysis presented here to include also the strange sector (kaons 
and eta) as well as to extend previous studies of transport coefficients by including pion chemical potentials along the 
lines presented here. 
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APPENDIX A: HOLOMORPHIC PATH INTEGRALS 

We review here some of the key aspects of the holomorphic path integral representation which are used in the main 
text. We will follow the discussion in (53], to which we refer for more details. 

We consider the space S of complex analytic functions of one complex variable and define the following scalar 
product: 

(f\g) =1%^- e-*> W)9{z) (Al) 

where the bar denotes complex conjugation (z and z are treated as independent variables), and the notation for the 
measure means 

oo 

dz dz f dx dv , . . 

^r~= / A2 

Zm J 7r 

— oo 

with z = x + iy. We also define the states (z\ in the dual space S* such that (z\f) = f(z), with |/) e S. Then, the 
set {/„}§°, with 

fn(z) = 4= (A3) 



constitutes an orthonormal basis for S with the inner product ( Al I. This implies in particular: 

dz' dz' 



2m 



e" z z e z z f(z') = f(z) (A4) 



Wc can also calculate the scalar product: 

oo oo 1 

z') = Uz)W) = £ -y(zzT = e** (A5) 



Z2 



n=0 n=0 



where we will denote the dual of (z\ by \z). Now, from the definition (Al I, the identity operator can be written as: 



1= l%^e---\z){z\ (A6) 
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Since the functions (A3 1 constitute an orthonormal basis, we can calculate the trace of an operator as follows: 



n=0 



dz dz 
2iri 



(A7) 



The prescription ^ defines a representation of the creation and annihilation operators on S. Therefore, 

/ i i _/\ & i i _/\ _/ 



{z\a)\z') =z(z\z') = zc z 



(Z\a\z 



(A8) 



For the purpose of obtaining a path integral, we need to know how to calculate the matrix elements (kernels) of the 
kind 0(z, z') = (z\0(a\ a)\z'), where O is an operator expressed in terms of creation and annihilation operators. If 
the operator is expressed in normal-order form (which we denote by the subscript N) i.e, arranging creation operators 
to the left and annihilation to the right, the kernels can be written in a particularly useful way, from (A8), as: 



(z\6 N \z') = N (z,d/dz)e zS ' = N (z,z')e z 
In particular, we will need the kernel corresponding to "time-evolution": 



(A9) 



U(z,z';t f -U) = (zle-^f-^lz') 



(A10) 



with £/n the normal-ordered hamiltonian of the system. 

For that purpose, as customarily, we divide the interval tf — U into n subintervals of infinitesimal length e and we 
will take the n — > oo limit in the end. Having in mind the application to thermal field theory, we will take complex 
times t S C where C is the contour starting at ti and ending at tf = ti — i/3 showed in FigJTJ 

Now, from (A9l, for an infinitesimal time interval: 



U(zi,z 2 ;s) ~ ^h^m)^ 



(AH) 



so that, inserting the identity operator (A6l n—1 times in (A10) one gets: 



/n— 1 , j_ f n—l 

n ^27r/ fc 6XP I ZlE ' + ^ [( Zk+1 ~ Zk ^ k + £ ^N(z k+ll z k ) 
k=l m I fc=l 



(A12) 



with z n — z. 

For the Hamiltonian ^ one has H = Hn + lu/2 and the previous integral can be explicitly calculated by using the 
standard formula [23] : 



n 

k=l 



dz k dz k 
2ni 



-zAz-\-uz-\-uz 



(det A)~ e 



l„uA~ 1 u 



(A13) 



which, taking the n —> oo limit, yields: 



U (z,z';t f - U) = exp (zz'e'^f-^ + S\j]\ 
where the subscript "0" distinguishes the particular case of the hamiltonian ^ and 



m = * j dt 

c 



-iu>(t f -t) 



z- ^— j(t)+z'- 



j(U-t) 



>2uj 



I2u 



dt dt' j(t)8(t-t')- 



o(t~t') 



2oj 



(A14) 



(A15) 



25 



APPENDIX B: FREE THERMAL PROPAGATORS AND PARTITION FUNCTION AT n ± 

In this appendix we review some important aspects regarding the canonical description of the free theory and the 
different representations for free propagators in Thermal Field Theory at nonzero chemical potential, paying special 
attention to the differences between the case of particle number chemical potential and that of exactly conserved 
charges such as the electric charge for complex scalar fields. 

Let us consider first the case of a free neutral scalar field 4>{x). In that case, one can evaluate the partition function 
and the propagator (two-point function) directly in the complete set \Ni, N 2 , . . .), corresponding to eigenstates of the 
hamiltonian operator with N\ particles in the state 1, N 2 particles in the state 2, and so on: 

(N U N 2 , . . . \N[,N 2 , ...)= 6 NlN , ■ S N2N , 

/ °° E \ 

Ho\N!,N2,...) = ( N 1 E 1 + ... + J2y) \Ni,N»,...) 

N\N 1 ,N 2 ,...) = (N 1 + N 2 + ...)\N 1 ,N 2 ,...) (Bl) 

with ^2 Ni = N. For non-interacting bosons of mass to, E; k = \J to 2 + \pi\ 2 , and the (infinite) term y\- Ej/2 is the 
vacuum energy. As customary, we consider first the system in a finite volume V — L 3 , which we will later take 
to infinity, so that spatial momenta are discretized as \pi\ = ^jf- with integers rii and energy levels are labeled by 
ft = (n x ,n y ,n z ). The free partition function reads then: 

oo —BEa/2 

n = n e e-^-M) e -^ = n ^v^, 

ft N=0 n 

where the condition /i < Eft must be satisfied for all n. Thus, in the V — > oo limit: 



logZ° g = -V 



d 3 p 
(2tt) 3 



PE p 



log 1 



(B2) 



where E 2 = \p\ 2 + to 2 . Therefore, in the following we must restrict to a chemical potential fi < to (below the 
Bose-Einstein condensation limit) to ensure the convergence of the previous expressions. 

In order to obtain the free particle propagator in the canonical formalism, defined as the two-point function: 



G(x) ^ (f$(x)$(O)) 0!fi = Z- 1 Tr{e-^-^>T $(x)$(0)} 



(B3) 



where T is the time-ordering operator, we expand the field as customarily in terms of creation and annihilation 
operators: 



^ilirri-x I L _|_ ^ — i2-Kn-x / L 



(B4) 



V ^ y/2E; 

n v 

with commutation relation 

[^,&U = V5 H ^ (B5) 
The free hamiltonian and the number operator are given in terms of creation and annihilation operators as: 

ft 

^ = E^4«™ (B6) 



Now, the reaZtime evolution of the field is given by <p{t,x) = e lHt <j)(x)e lHt with igl. We will calculate the trace 
in ( B3 1 using 
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so that we get for the free propagator, after taking the V — > oo limit: 

G(x) = 9{t)G>{x) + 9(-t)G<(x) (B8) 

with: 

G>«)(x)= | ^ e * 2 G>«)(t,p) (B9) 

G>(t,p) =^ [e-^ ( (l + n(£ p - A 4 )) + e iEpt n(E p - p)] (BIO) 

G<(t,p) =^ [e^*(l + n(£ p - fi)) + e- iEpt n(E p - „)] (Bll) 

Note that the above propagators are obtained from the p = ones by the following replacement in the distribution 
function: 

n(x) -> ra p (x) ee ^— 1 (B12) 

with: 

^ = 13 i 1 ~ i~) (B13) 

Therefore, we have for instance h p {E p ) = n(E p — p) and the h p function satisfies: 

1 + h p (x) + h p (-x) = (B14) 
Thus, the free propagator satisfies the following KMS-like periodicity condition in the mixed representation: 

G>(t,p) = G<(t + z/3 p ,p) (B15) 
and in Fourier space we can write a spectral representation: 

G > (p ,p) = [1 + %(p )] p(po,p) 

G<(p ,p) =e-^°G>(p ,p) = h p ( Po )p(p ,p) (B16) 

where 

p(po,p) = 2™ g n(po)<5 [(p ) 2 - E 2 p ] (B17) 

is the free spectral function, which is independent of temperature and chemical potential. 
Now, using: 

-oo 2tt k + te 



with e — > , we can write for the propagator in (B8) in momentum space 
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FIG. 10: Contours used to derive the Matsubara representation of the free imaginary-time propagators. The black dots on 
the imaginary axis denote the Matsubara frequencies and 7? — *• oo. 



G(p ,p) 



Vt-El 



•27r*(pg - E 2 p )n(\p \ - n) 



(B19) 



Note that we have used n p (E p )S(pQ — E p ) = h p (\p \)5(pQ — E p ) = n(\p \ — p)S(pq — E p ) and we have chosen the 
"|p |-prescription" which, as explained in the main text, guarantees the decoupling of the imaginary-leg contribution 
to real-time Green functions. 



The free propagators in (BIO I-(|B11 1 can be extended to imaginary times t — —it corresponding time differences 



along the imaginary-time leg C4 in Fig l| Thus, we define At(t,p) — G > (—ir,p) for r > and At(t,p) = G < (—iT,p) 
for t < 0. Now, if we try to construct a Matsubara frequency representation in this case, we have, from the mixed 
representation (BIO I and using (B12 l-(B13 1 : 



A T (r > 0,p) 



1 

2~7ri 



C1UC2 e 



/3 p z 



1 z 2 -E* p 



a<T<p p 1 



El 



(B20) 



where the contours are showed in Fig 10 the black dots on the imaginary axis being the "modified" Matsubara 
frequencies cu n = 2nn/(3 p . A very important point here is that the last step in the above equation is only valid for 
r e [0,/5 p ], otherwise the integrals along the circular arcs with R — > 00 do not vanish. Thus, the Matsubara Fourier 
representation is only valid in that interval, which is smaller th an [0, 0\. 

Carrying out the same procedure with At(t < 0,p) using (Bill leads to th e same "modified" Matsubara rep- 
resentation for r G [— j8p,0]. In fact, we see that the KMS-like condition ( B15 1 translates into the imaginary-time 
propagator as: 



A T (r + /3 p ,p) = A T (r,p) 



(B21) 



so that this propagator does not satisfy the usual equilibrium KMS condition At(t + (3,p) = Ar(r,p). 

At this point, it is instructive to compare the above free propagators for chemical nonequilibrium with those obtained 
when an exact conserved charge is present. As we are going to see, there are crucial differences between the two cases. 
For definiteness, we consider the electric charge for the ca se of a com plex scalar field and denote the corresponding 
chemical potential by /iq. In that case, the counterparts of (BIO I-(B11 1 for the free propagator Gq(x) = (T eft (x) (f)(0)) 
are [TTj : 



G%if,p) =2£- [e- iE *'(l + n(E p - MQ )) + ^n(E p + (j.q)] 



2E n 



[e iE p*(l + n(E p + [xq)) + e- zE ^n(E p - n Q 



(B22) 
(B23) 



Note that, unlike our previous case in (BlO)-(Bll), the chemical potential enters now with opposite sign for the 
positive and negative frequencies, which comes from the opposite charge of particles and antiparticles, necessary to 
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maintain the chemical equilibrium imposed by charge conservation. Due to this sign difference, the above propagator 
satisfies now the following condition: 



GZ(t,p) = e-^G<(t+if3,p) 



(B24) 



i.e, in this case the KMS symmetry realizes simply as a modification of the fig = KMS boundary condition by a 
constant /XQ-dependent factor, which is a result completely different from the previous case, c.f., eq.(B15l, where the 
loss of KMS involves the p-dependent /3„, which cannot be rewritten as a multiplicative factor: 



G>(t,p) = 



1 

2E„ 



-^ e -i Ep (t+i0) n ^ Ep _ ^ + e 0^ e iE p (t+i0)^ + n ( Ep _ ^1 ^ e -fi^G<{t + i0,p) (B25) 



Another way to arrive to the same conclusion is to derive directly the periodicity relation from the thermal averages. 
In the charged scalar case, it is crucial to use that the field is a charge eigenstate, i.e. [Q, <j)] = —(f), [Q, <ft] = ft |17j . 
This, together with charge conservation [Q,H] — leads to: 



Tr 



(B26) 



However, in the case of a real field and the number operator, even though [N, H) = in the free case, [N, <fi] ±cf> 
which prevents the previous relation to hold. 

Defining now the imaginary-time propagators as above, we get the same factor in the /xq case: 



A?(r + Ap)=e-^A«(r ( p) 



(B27) 



In fact, it is not difficult to see that in the /xq case, this simple form of KMS symmetry still allows for a well-defined 
Matsubara IT frequency representation: 



2ni 



0<r</3 1 



C'UC' e 



/3(z+MQ) 



1 Z 2 -£2 



(B28) 



where C' 12 correspond to the contours in Fig 10 but with the vertical line displaced to z — (hq < m) and the 

dots in that line being now the standard Matsubara frequencies uj n — 2-Kn/(3. Therefore, in this case the ordinary IT 
formalism is recovered for r 6 [— [3, [3] simply by changing in the Feynman rules uj n — > uj n + ifiQ. 

Most of the results showed in the main text can be written in terms of the above thermal propagators evaluated at 
the origin in position space and functions related to them. From (B9 1-( Bl 1 1 we have (for \x < m) at r = t = x = 0: 



G>(0) = G<(0) = A T (0) = G(0) = G(0) 



T=/x=0 



gi(m,T, n) 



(B29) 



where the T — /i = contribution is ultraviolet divergent. In dimensional regularization it is given by: 



G(0) 



T=n=0 



jD-1 



P 1 



r[i 



(2tt) d ~ 1 2E p (4tt) d / 2 



(B30) 



while the T, ^-dependent contribution gi is finite. We are following the same notation as in |31j so that g\ is the /x ^ 
extension of their function gi(T), to which it reduces for /i = 0. We have: 



gi(m,T,fi) = 



2tt 2 



dp 



E p e^p-ri - 1 



(B31) 



Note that in dimensional regularization one has, as in the /x = case: 
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[d 2 -V 2 ] A T (r,f) 
OG(t,x) 

and df,G(0) = <9 M A T (0) = 0. 

Let us also define, following again the notation in [3"T] : 



m 2 A T (0) 
~m 2 G(0) 



(B32) 



g (m,T,n) 



rp />00 



7T 



dp p log 



1 - e 



-/3(B P -M) 



(B33) 



so that, taking into account that dE p /dm 2 = l/(2E p ), we can write the free partition function (B2) separating its 
divergent contribution in dimensional regularization as: 



log Z% = 



(iV 



(4tt) d / 2 



+ g (m,T, (j,) 



(B34) 



Finally, note that the functions go and gi satisfy a similar relation as in 2 



<7i(m,T, y) 



dm 2 



go(m,T, fj.) 



(B35) 
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